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ABSTRACT 

We describe Hubble Space Telescope (HST) imaging of 10 of the 20 ESO Distant Cluster Survey 
(EDisCS) fields. Each ~40 square arcminute field was imaged in the F814W filter with the Advanced 
Camera for Surveys Wide Field Camera. Based on these data, we present visual morphological 
classifications for the ~920 sources per field that are brighter than I au t = 23 mag. We use these 
classifications to quantify the morphological content of 10 intermediate-redshift (0.5 < z < 0.8) galaxy 
clusters within the HST survey region. The EDisCS results, combined with previously published data 
from seven higher redshift clusters, show no statistically significant evidence for evolution in the 
mean fractions of elliptical, SO, and late-type (Sp+Irr) galaxies in clusters over the redshift range 
0.5 < z < 1.2. In contrast, existing studies of lower redshift clusters have revealed a factor of ^2 
increase in the typical SO fraction between z = 0.4 and z = 0, accompanied by a commensurate 
decrease in the Sp+Irr fraction and no evolution in the elliptical fraction. The EDisCS clusters 
demonstrate that cluster morphological fractions plateau beyond z w 0.4. They also exhibit a mild 
correlation between morphological content and cluster velocity dispersion, highlighting the importance 
of careful sample selection in evaluating evolution. We discuss these findings in the context of a 
recently proposed scenario in which the fractions of passive (E,S0) and star-forming (Sp,Irr) galaxies 
are determined primarily by the growth history of clusters. 

Subject headings: galaxies: clusters: general — galaxies: formation — galaxies: evolution 
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Although the morphology-density relation is observa- 
tionally well-established, it is still unclear why the inci- 
dence of early type (E+S0) galaxies is higher in over- 
dense regions of the universe. The options are often 
couched in terms of "nature" versus "nurture". These 
scenarios describe different paths to qualitatively similar 
relations between morphology and environment, and are 
therefore difficult to disentangle. In the nature scenario, 
galaxies that end up in high density environments at low 
redshift are more likely to have experienced initial con- 
ditions leading to an early-type morphology upon forma- 
tion. In the "nurture" scenario, all galaxies, regardless of 
their local densities at low redshift, have identical proba- 
bilities of having formed as early types. The morphology- 
density relation is then the result of subsequent morpho- 
logical alterations as galaxies enter increasingly higher 
density environments as structure grows. Many mech- 
anisms that could produce the required morphological 
transformations have been suggeste d, including merg- 
ers and galaxy - galaxy interactions ( j Tbomre fc Toomrg 
[19721: [Tck3 [l985l iLaverv fc HenrvHl988tlMihodl2004D . ha- 



rass ment d Richstonc 1976; Moore ct al. 1998), gas strip- 
ping (|Gunn fc Gotdll97llAbadi et al 1999; Quilis et al] 
2000), strangulation (|Larson et al.l |1980T iBekki et al.l 



2002), and cluster tidal forces (|Bvrd fc Valtonenl ll990). 



The observed evolution of the morphology-density re- 
lation provides an important constraint on models of its 
origin. Observations of the high-density regions of clus- 
ter cores show that while the elliptical fraction has not 
evolved, the SO (Sp+Irr) fraction has grow n (diminished) 
by a factor of ~2 over the last 5 Gyr (Dressie r et al.1 
Il997t iFasano et al.l 2000). However, over the same time 
interval, no evo l ution is observed at lower densities 
(|Treu et all 120031: ISmith et al.ll2005l) . The morphology- 
density relation at larger lookback times (z > 0.5) has 
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just begun to be explored. ISmith et al.1 (|2005f) used Hub- 
ble Space Telescope (HST) Wide-Field Planetary Cam- 
era 2 (WFPC2) data for six clusters at 0.75 < z < 1.25 
to estimate the early-type fraction as a function of local 
galaxy surface density. They extended their measure- 
ments to low densities (~1 Mpc -1 ) by including field 
galaxies in HST WFPC2 images of the cluster CL0024 
at z = 0.395. Com paring th eir results to similar studies 
at lower redshift s (jPresslerl Il980bt iDressler et all 119971 ; 
iTreu et ail |2003[ ) , they found that the early- type frac- 
tion has increased steadily since z ~ 1 in the densest 
regions (^1000 Mpc -2 ), has increased only since z ~ 0.5 
in moderate density regions (~100 Mpc -2 ), and has re- 
mained consta nt since z ~ 1 in the lo west density regions 
(<10 Mpc" 2 ). iPostman et~aTl (|2005f) find compatible re- 
sults in their study of the morphology-density relation 
using Advanced Camera for Surveys (ACS) HS T imag- 
ing of seven z ~ 1 clusters. ISmith et al.l (|2005h suggest 
a nurture scenario to explain the observed evolution in 
the morphology-density relation. If structure growth is 
hierarchical, the densest regions at any redshift collapsed 
the earliest. Thus, if environmental processes can modify 
a universal initial mix of morphologies through a trans- 
formation from late to early types, these processes would 
have had a longer time to operate in increasingly dense 
regions. In addition, the efficiency of such transforma- 
tions could be density-dependent. 

To discover whether environmental processes, rather 
than initial conditions, drive the differential evolution of 
the morphology-density relation, we must first demon- 
strate that nature scenarios cannot reproduce the obser- 
vations. If we find that environmental factors are indeed 
important, it is critical to identify the responsible mech- 
anisms. The current data are insufficient to discrimi- 
nate among the alternative models. For this, we require 
a full mapping between galaxy properties and environ- 
ment. Towards this goal, we present the morphological 
fractions in ten rich clusters at 0.5 < z < 0.8 . drawn from 
the E SQ Distant Cluster Survey 16 (EDisCS: (White et all 
2005). Because strong evolution has been observed in 
high-density regions, clusters are a critical environment 
to probe. The redshift ran ge of o u r study bridges a ga p 
bet ween those of iDresslerl (|1980eJ ). iFasano et ail (2000), 



and IDressler et ail (|1997) at low redsh if t and those of 



Postman et~afl ([2005) at high 



ISmith et alJ ([20051 ) and 
redshift. The EDisCS clusters represent a significant 
increase in the number of well-studied clusters at high 
redshift. Because the mean age of clusters decreases 
with increasing redshift, the scatter in cluster properties 
may also increase with redshift. Thus, large samples are 
necessary to disentangle evolution from cluster-to-cluster 
variations. Finally, we measure morphologies from high- 
quality HST ACS images, allowing us to distinguish SOs 
from other early types. As described above, the SO pop- 
ulation appears to have increased at the expense of the 
Sp+Irr population since z = 0.5. Thus, tracing the SO 
population to high redshifts may place important con- 
straints on their formation. 

This paper is organized as follows. In iJ5]we describe 
the EDisCS sample. In g3]we describe the HST ACS data 
used to measure morphologies. Our classification proce- 
dure and Hubble types for all galaxies with J au t < 23 

16 http: //www.mpa-garching.mpg. de/galf orm/ ediscs/index . shtml 



mag are presented in []4j In $5] we describe our method 
of quantifying the morphological content in each clus- 
ter. Our results are presented in SJ5J and our conclusions 
regarding them can be found in 

Results will be presented for two sets of cosmological 
parameters, to allow a direct comparison with previous 
work (fi = 1, A = 0, H = 50 km s -1 Mpc" 1 ) and 
to provide the corresponding estimate for the currently 
standard cosmology (fio — 0.3, A = 0.7, Ho = 70 km s -1 
Mpc -1 ). In the remainder of this paper we will refer to 
these as the classic cosmology and the WMAP cosmol- 
ogy, respectively. 

2. THE ESO DISTANT CLUSTER SURVEY (EDISCS) 

The ten clusters included in this study are a subset of 
the ESO Distant Cluster Survey, an extensive program 
to obtain spectroscopy and multiwavelength photometry 
for 20 galaxy clusters at 0.4 < z < 1. Details of the goals 
an d data sets as s ociate d with this survey are documented 
bv lWhite etafl (|2005h . but we briefly summarize below. 

The clusters chosen to be part of EDisCS were 
drawn from the Las Campanas Distant Cluster Survey 
(LCDCS) catalog, which consists of optically-selected 
clusters identified from fluctuations in the extragalactic 
background light (iDalc antonl 119961: IZaritskv et~aJU Il997t 
iGonzalez et al.ll200lD . Theredshifts of the LCDCS clus- 
ter candidates were initially estimated from the apparent 
magnitude of the brightest cluster galaxy (BCG). The 
candidates were divided into an intermediate-redshift 
group (z ~ 0.5) and a high-redshift group (z ~ 0.8). 
Shallow two-color imaging was obtained at the Very 
Large Telescope (VLT) for the 30 brightest candidates 
in each group. A sample of 20 EDisCS clusters was 
selected based on the presence of a re d sequence and 
a BC G of the appropriate magnitude ([Gonzalez et al.l 
2001). These twenty clusters were then imaged more 
deeply at the VLT in BVI (45 minutes) for the z — 0.5 
candidates and in VRI (2 hours) for the z ~ 0.8 can- 
didates. These data, covering a 6.5' x 6 .5' region in 
the fie ld of each cluster, are presented in IWhite et al.l 
(2005), an d form the basis of t he weak-shear analysis pre- 
sented by [Cloweitan dHH) . Additional near-infrared 
imaging was obtained with SOFI at the New Technology 
Telescope (NTT) in K s for the z ~ 0.5 candidates and 
in JK S for the z ~ 0.8 candidates (Aragon-Salamanca 
et al., in preparation). The near-infrared imaging covers 
approximately 6' x 4.2' in the fields of the intermediate- 
redshift clusters and 5.4' x 4.2' in the fields of the high- 
redshift clusters. The optical and near-infrared imaging 
were used to measure photometric redshifts, as described 
in Pcllo et al. (in preparation). 

An initial phase of spectroscopy consisted of relatively 
short exposures of a single slit mask per cluster. One 
cluster with an estimated redshift of z ~ 0.8 was re- 
vealed as a superposition of weak groups, and was re- 
jected from the sample. These spectra revealed that 
the remaining 19 clusters have redshifts in the range 
0-4 % z < 1. Deeper spectroscopic exposures (~2 hours 
each for galaxies with 18.6 < I(r = 1") < 22 mag in 
intermediate-redshift clusters; ~4 hours each for galax- 
ies with 19.5 < I(r = 1") < 23 mag in high-redshift 
clusters) of 3-5 masks per field were then taken, and 
are of sufficient quality to obtain information about the 
stellar populations and internal dynamics of the target 
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ga laxies. The sp e ctrosc opic observations are described 
in lHallidav et al.l (|2004l ) and Milvang-Jensen et al. (in 
preparation) . 

The ground-based EDisCS imaging and spectroscopy 
go a long way towards characterizing both the clusters 
themselves and the galaxies within them. At the high 
redshifts of the EDisCS sample, however, only HST can 
provide the spatial resolution necessary to provide robust 
morphologies. Motivated by the issues discussed in SJTJ 
we obtained HST imaging for the ten highest redshift 
EDisCS clusters, listed in Table [1] along with their basic 
physical parameters. Cluster 9 (cll232-1250) was drawn 
from the EDisCS intermediate-redshift sample, and the 
remainder were drawn from the high-redshift sample. 

3. HST ACS DATA 

The HST observations were designed to coincide as 
closely as possible with the coverage of the ground-based 
optical imaging and spectroscopy, within guide star con- 
straints. The ground-based optical data cover a 6.5' x 
6.5' region around each cluster, with the cluster center 
displaced by 1' from the center of the region. For ref- 
erence, the ACS Wide Field Camera has a field of view 
of roughly 3.5' x 3.5'. Balancing scientific motives for 
going deep over the entire field against a limited num- 
ber of available orbits, we tiled each 6.5' x 6.5' field 
in four pointings, with one additional deep pointing on 
the cluster center (taken as the location of the BCG). 
The resulting exposure time per pixel is 2,040 seconds, 
except for the central 3.5' x 3.5', which has an expo- 
sure time per pixel of 10,200 seconds. The deep central 
pointing probes to lower surface brightness, fainter mag- 
nitudes, and larger galactic radii in the region of the 
cluster containing the most galaxies. The centers of the 
structures ultimately chosen for spectroscopic followup 
in the clll03-1245b and cll227-1138 fields are somewhat 
offset from those anticipated at the time the HST ob- 
servations were taken. All exposures were taken under 
LOW SKY conditions to maximize our surface brightness 
sensitivity. 

The ACS calibration pipeline, CALACS v4.3 (6- 
June-2003) debiased, dark-subtracted, and flat-fielded 
our ACS images "on the fly" when they were requested 
from the HST archives. Known bad pixels and satu- 
rated data were also flagged in accompanying data qual- 
ity (DQ) arrays. Approximate World Coordinate Sys- 
tem (WCS) headers were provided. The resulting images 
were returned by the archive with the FLT suffix. In ad- 
dition, images that were cosmic-ray split were combined 
and returned with a CRJ suffix. 

To produce a mosaic, we required precise offsets be- 
tween the 32 FLT images retrieved for each cluster. The 
approximate WCS headers provided were insufficient for 
this purpose. We found that the offsets between cosmic- 
ray split images were negligible. We therefore computed 
the necessary shifts between pointings using the higher 
signal-to-noise CRJ images. The shifts could not be com- 
puted using the retrieved CRJ files because they had 
not been undistorted. We therefore drizzled each CRJ 
image separately, using the approximate WCS headers 
provided. Because each undistorted CRJ image overlaps 
the central pointing, one of the images centered on the 
cluster center was chosen as the reference image. Shifts 
between the reference image and each undistorted CRJ 



image were then computed using cross-correlation. 

Combination of the FLT images using the resulting 
shifts was accomplished using MultiDrizzle, a Python 
code written by Anton Koekemoer to run under PyRAF, 
the Python-based interface to the Image Reduction and 
Analysis Facility (IRAF) . MultiDrizzle automatically re- 
moves cosmic rays and combines dithered images using 
PyDrizzle, which has been developed by the Science Soft- 
ware Branch at the Space Telescope Science Institute. 
For each of the 32 FLT images per cluster, MultiDrizzle 
included negative bad pixels in the data quality array, 
subtracted the sky, and separately drizzled and undis- 
torted each image. Next, it created a median image from 
these separately drizzled images using shifts computed 
from the headers along with the user-supplied refinement 
shifts described above. This median image, relatively 
free of cosmic rays, was compared to the input images 
to identify cosmic rays. The median image was then re- 
distorted to create cosmic ray masks. These masks were 
used in the final image combination step using drizzle 
and the lanczosS kernel, which provided optimal noise 
properties. 

4. GALAXY MORPHOLOGIES 

We visually classified all galaxies bri ghter than J au to = 
23 m ag. Here I a uto is the SExtractor (|Bertin fc Arnouti 
119961 ) AUTO magnitude measured on the I-band VLT 
images, and is an estimate of the total magnitude of a 
galaxy, in the Vega system. This magnitude does not 
include an aperture correction, nor has it been corrected 
for Galactic extinction. The limit was set both to ensure 
robust classifications and to provide a tractable sample 
of ~9200 galaxies. 

Our classifications arc most useful if they conform 
to systems adopted by previous studies. For this rea- 
son, each classifier trained on the HST WFPC2 images 
and visual morphological catalogs of the 0.3 < z < 0.5 
M ORPHS clusters, using the same procedure described 
in lSmail et al.1 (|1997f ). For uniformity, each classifier used 
the same IRAF script to examine and classify EDisCS 
galaxies. This script displays two side-by-side versions of 
a 200 x 200 square pixel cutout centered on each galaxy 
meeting the magnitude limit described above. One ver- 
sion is on a log scale between -0.1 - 2 DN/s, while the 
other is on a log scale between -0.1 - 25 DN/s. Together, 
these displays allowed classifiers to inspect the galaxies 
from their high surface brightness cores to their low sur- 
face brightness outer features. In general, we found that 
the depth and quality of the EDisCS ACS data were 
similar to or better than the MORPHS WFPC2 data for 
lower redshift clusters. 

Classifications were performed by 5 of the authors 
(A AS, JJD, VD, PJ, BP). Each classified the galaxies 
in 3-6 clusters. For a given galaxy, the final Hubble type 
was based on the classifications of two or more of the 
authors. First, the Revised Hubble Type was translated 
into a T-type according to the scheme presented in Table 
[21 Classifications appended by one question mark were 
given half weight; those appended by two question marks 
were given one-quarter weight. If a classifier specified 
two types separated by a slash, the first was given three- 
quarters weight and the second was given one-quarter 
weight. All T types submitted for a given galaxy were 
then ranked by their frequency. The final Hubble Type 
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was then chosen from the highest ranked T types. In 
general, the highest ranked T type was chosen. If there 
were two equally-ranked classifications, one was chosen 
randomly If there were three or more equally-ranked 
classifications, the differences in T types were computed 
between them. If there were no gaps smaller than three 
T types, the median value was chosen. If there was one 
gap smaller than three T types, one of the two similar 
classifications was chosen at random. If there was more 
than one gap smaller than three T types, the median 
value of all classifications was chosen. 

The main results of this work are sensitive to how accu- 
rately galaxies can be placed into the broad categories of 
elliptical (E), SO, and late (Sp+Irr). Although accuracy 
in this case is difficult to quantify, we can test how con- 
sistently galaxies are placed in the same bin by different 
classifiers. Using ^900 galaxies down to / au to = 23 mag 
in cll216-1201, all of which were assigned a Hubble type 
by all classifiers, we computed the raw fractions for each 
classifer, uncorrected for the presence of foreground and 
background galaxies. The root mean variance for any 
morphological fraction is < 0.10, which is comparable to 
the error computed from Poisson statistics. 

There is some controversy surrounding the ability of 
morphologists to distinguish between E and SO galaxies, 
especially at high redshifts. For this reason, some pre- 
vious investigators chose to lump E and SO galaxies to- 
gether into an early- type class. For easy comparison with 
these works, in the following we have plotted and tabu- 
lated our results for early types as well as for E and SO 
galaxies separately. However, because the SO population 
is strongly evolving in clusters at z < 0.5, it is important 
to track this population at higher redshifts. Comparing 
the EDisCS ACS data to the MORPHS WFPC2 data 
gives the impression that we can distinguish SOs from el- 
lipticals with the same accuracy for the EDisCS sample 
as was obtained with the MORPHS sample. To test this 
impression, we have performed a quantitative analysis of 
the light distributions of E and SO galaxies within cll216- 
1201, which is our highest-redshift cluster and which has 
been visually examined by all classifiers. Specifically, 
we used the ELLIPSE task in the Space Telescope Sci- 
ence Data Analysis System (STSDAS) package of IRAF 
to measure the surface brightness of a series of isopho- 
tal ellipses in each galaxy as a function of semi-major 
axis. For each step, the semi-major axis is increased by 
a factor of 1.1. The ELLIPSE task was allowed to vary 
the ellipticity and position angle of the ellipses within a 
given galaxy. For each surface brightness profile, we fit 
a bulge+disk model of the form 

ibulge(r) = I e exp (-7.67[(r/r e ) 1/4 - 1]) (1) 

and 

Idisk(r) = Io exp (-r/r<i), (2) 

where the total surface brightness is given by /bulge + 
/disk- We considered three parameters for which E and 
SO galaxies should display different distributions: 1) the 
width of the ellipticity distribution in a given galaxy, 
R, defined as the interquartile range of all values of the 
ellipticity for a given galaxy; 2) the bulge-to-total frac- 
tion B/T; and 3) x 2 goodness-of-fit from fitting the r 1 / 4 



law (i.e. only bulge, no disk). The differential and frac- 
tional cumulative distribution functions for these three 
parameters are plotted in Figure [TJ At a statistically 
significant level, SOs display larger ellipticity distribu- 
tion widths than ellipticals. The presence of a disk in a 
galaxy will result in a wider range of ellipticities, since 
the ellipticity of the galaxy changes as the disk begins 
to dominate the surface brightness profile. Likewise, the 
B/T fractions for visually-classified ellipticals appear to 
be skewed to higher values compared to SOs, as expected. 
This is confirmed by a two-sided Kolmogorov-Smirnov 
test. Finally, a higher fraction of ellipticals have small 
values of x 2 , indicating a better fit to a pure bulge model. 
However, the differences in the x 2 distribution are not 
significant. The overall results of these tests are that 
visually-classified SOs have quantitatively different sur- 
face brightness profiles than visually-classified ellipticals. 
Although in some individual cases it is difficult for visual 
classifiers to differentiate between Es and SOs, it is clear 
that, in a statistical sense, objects visually classified as Es 
and SOs form two distinct populations with objectively 
measurable physical differences. More sophisticated two- 
dimensional profile fits will be presented in Simard et al. , 
in preparation. 

Morphological catalogs for the ten highest-redshift 
EDisCS clusters are available in the on-line version of 
this paper, with columns described in Table O 

5. ANALYSIS 

The goal of this paper is to quantify the morpholog- 
ical content of the EDisCS clusters. In particular, we 
will discuss the overall fractions of E, SO, and late-type 
galaxies in our clusters and compare them to the frac- 
tions f ound in clusters spann ing a range of redshifts. Fol- 
lowing Dressie r" et al.l (|1997l ). we do not include galaxies 
that were unclassifiable (types -6 and 66 in Table [5]) in 
our analysis. To facilitate a fair comparison with other 
samples, we compute the morphological fractions as con- 
sistently as possible with previous work. In the following 
subsections, we describe the key elements of our analysis. 

5.1. Magnitude Range 

Early types preferentially occupy the bright end of the 
galaxy luminosity function, while late types dominate 
the faint end (e.g. iBlanton et aLll2001t iGoto et al.ll2002t 
IZucca et al.l 120061) . Thus, morphological fractions de- 
pend upon the range of absolute magnitudes sampled. 
Morphological fractions in the MORPHS clusters were 
determined using apparent magnitude cuts in Isi4 de- 
signed to correspo nd to a V-band abso lute magnitude of 
My = —20 mag ([Dressier et al.l fl997l) . but in actuality 
corres ponding to My = — 19 mag d ue to a transc r iption 
error dFasano et al.l 2000). However. IFasano et alj (|2000h 
re-analyzed the MORPHS data with the intended limit- 
ing absolute magnitude and additionally analyzed nine 
clusters with 0.1 < z < 0.25 in the same manner, pro- 
viding an ideal comparison sample for the EDisCS clus- 
ters. We therefore adopt an absolute magnitude limit 
of My = —20 mag. For a given cluster, the limiting 
apparent magnitude in the /-band is then given by 

/auto,iim = Mv+51og 10 (d L , pc )-5-(Mv-Mi) + fci, (3) 
where d^ tPC is the luminosity distance of the cluster in 
parsecs (calculated in either the classic or WMAP cos- 
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mology, as indicated), My — Mj is the rest- frame color, 
and ki is the J-band k-correction. We adopt the rest- 
frame color and Cousins /-ban d k-correc t ions fo r an el- 
liptical galaxy, as presented in iPoggianll (|1997f ) . Equa- 
tion [3] results in values of / a uto,iim ranging from 21.3 to 
23.3 mag, the faintest values being slightly fainter than 
the limit of our visual classifications. Given the errors in 
the k-corrections relative to these differences, we adopt 
/ au to,iim = 23 mag in these instances. This occurs for 
cll216-1201 when using the WMAP cosmology and for 
cll054-1245, cll216-1201, and cll354-1230 when using 
the classic cosmology. 

5.2. Aperture 

Clusters exhibit a radial gradient in their mor- 
phological fractions; the centers of clusters con- 
tain a lar ger fraction of el l iptica l s than the out- 
skirts (e.g. iMelnick fc Sargentl Il977t iGoto et alj 12004 
iThomas fc Katgertll2006ir Thus, in comparing the over- 
all morphological fractions across clusters, one must 
choose a consistent aperture. The MORPHS project has 
set the standard for the aperture within which to com- 
pute morphological fractions. Most studies of the mor- 
phological fractions of galaxy clusters at other redshifts 
have used a similarly-sized aperture. So that we may 
assess the level of evolution between the morphological 
studies presented in the literature and those conducted 
with EDisCS data, we adhere to this precedent and adopt 
a circular aperture of radius 600 kpc (classic cosmology) . 

Since the radial density profiles of clusters vary, the 
average galaxy density within a fixed metric aperture 
will also vary. We therefore also calculated the mor- 
phological fractions within an aperture that scales with 
-R20O1 the radius within which the average mass density is 
equal to 200 times the critical density. The d erived i^no 
values were computed using Equation 8 in iFinn et al.l 
(2005), and are listed in Table Q] for both the classic and 
WMAP cosmologies. Unfortunately, our imaging data is 
complete out to i? 2 oo for only cll040-1155, cll054-1146, 
cll054-1245, and cll354-1230. We therefore use a radius 
of 0.6 x i?200j which keeps the analysis area within the 
imaging region for all clusters except for cll227-1138 and 
cll232-1250. The fraction of the analysis region which is 
not included in these two clusters is very small, and the 
resulting effect on the morphological fractions is likely 
minimal. 

5.3. Background Subtraction 

Galaxies within the magnitude range and distance 
from the cluster center described in i j5.ll and ij5.2l lie at 
a variety of redshifts. We wish to determine the pro- 
portions of E, SO, and late-type galaxies of cluster mem- 
bers only. Because the morphological mix of field galax- 
ies differs substantially from that of cluster populations, 
morphological fractions computed without regard to field 
contamination will underestimate the fraction of early 
types and overestimate the fraction of late types in clus- 
ters. We use four methods, described in detail below, to 
account for the presence of cluster non-members, check- 
ing for consistency among the resulting morphological 
fractions. 

5.3.1. Spectroscopic Redshifts 



As described in Sj2j the EDisCS program includes an 
extensive spectroscopic survey. For each cluster, red- 
shifts were obtained for 6-86 galaxies which both meet 
our apparent magnitude limit and lie within 0.6 x -R200 
(WMAP cosmology). We use this spectroscopic sample 
to constrain the morphological fractions in each cluster 
in two ways. First, we calculate hard limits on the mor- 
phological fractions using both the spectroscopic sample 
and the purely photometric sample, consisting of galax- 
ies that were either not targeted for spectroscopy, or for 
which spectroscopy failed to yield a redshift. This calcu- 
lation does not require that the spectroscopic sample be 
complete. Second, we estimate the morphological frac- 
tions using only the spectroscopic sample after applying 
small completeness corrections. 

The upper and lower hard limits on the fraction of 
galaxies of type i are given by: 

and 

where the subscripts s and p refer to the spectroscopic 
sample and the photometric sample, respectively; and 
the superscripts m and u indicate cluster members and 
galaxies at unknown (spectroscopic) redshift, respec- 
tively. Each quantity is a function of galaxy type i. The 
fraction of galaxies of type i is minimized when all of the 
other types that do not have redshifts are members and 
all of the galaxies of type i that do not have redshifts are 
nonmembers. The fraction of galaxies of type i is max- 
imized when all of the galaxies of other types without 
redshifts are nonmembers and all of the galaxies of type 
i without redshifts are members. The results of applying 
equations H] and [5] are shown as the grey shaded regions 
in Figure^ 

To directly estimate the morphological fractions from 
the spectroscopic data, we must correct for incomplete- 
ness in the s pectroscopic sample. A s described in the 
Appendix of iPoggianti et ail J2006), the completeness 
of the EDisCS spectroscopy is a weak function of both 
apparent /-band magnitude and clustercentric distance. 
Using the cluster-by-cluster magnitude and geometric 
completeness weights W mag and W gco calculated in that 
work, we estimate the spectroscopic morphological frac- 
tions / S poc(i) as: 

aw*) = E ^magor^gcoor 1 (6) 

i=i 

A spoc (tot) = J2 wwor^geoO/r i . (7) 

/ S pec(i) = N spcc (i) / N spcc (tot) (8) 

The probability of measuring a given morphological 
fraction is the product of the Poisson probability of mea- 
suring the total number of cluster members times the 
binomial probability of measuring the observed number 
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of a given morphology. We use the approxim ation pre- 
sented in Equations 21 and 26 of lGehrela (|1986h to derive 
l-cr error estimates based on this premise, using A spe c(«) 
and N spec (tot) as inputs. The results of applying Equa- 
tion [8] are shown as hollow squares in Figure [2l 

As described in the faint magnitude limit of the 
EDisCS spectroscopic survey is I(r = 1") = 22 mag for 
the intermediate-redshift clusters and I(r = 1") = 23 
mag for the high-redshift clusters. In general, I(r = 1") 
is fainter than 7 au t . For six of our clusters, all of the 
galaxies meeting the I a uto magnitude and aperture re- 
quirements of our analysis have values of I(r = 1") 
brighter than the spectroscopic survey limit. However, 
for cll054-1245, cll216-1201, cll232-1250, and cll354- 
1230, 10-20% are fainter than the spectroscopic flux 
limit. Most (>80%) of these are late-type galaxies. Thus, 
it is possible that our direct spectroscopic estimates are 
biased slightly towards low late-type fractions. In Fig- 
ure [2j we compare the morphological fractions derived 
from different methods, and find that for cll054-1245, 
cll216-1201, and cll232-1250, the spectroscopic method 
does yield lower late-type fractions than the photometric 
redshift or statistical background subtraction methods. 
However, all methods produce late-type fractions that 
are consistent with one another, within the errors. For 
cll354-1230, the spectroscopic method results in a late- 
type fraction between the photometric redshift method 
and the statistical background method. 

5.3.2. Photometric Redshifts 

Our optical and near-infrared imaging allows the 
derivation of photometric redshifts, as described in Pello 
et al. (in preparation). Briefly, two estimates of the 
redshift probability distribution (P(z)) were computed 
for each galaxy. Two independent co des were employed, 
one de scribed in lRudnick et al.l (120011) andlRudnick et al.l 
(|2003| ) and Hyperz, described in lBolzonella et alJ (|2000). 
An estimate of P(z) can be integrated over a suitable 
interval Az (in this case ±0.1) around the cluster red- 
shift to obtain the probability P c i us t that the galaxy is a 
member of the cluster. We used our large spectroscopic 
sample to determine /"thresh, the minimum value of P c iust 
required for a galaxy to be considered a cluster member. 
Membership information derived from each of our two 
estimates of P(z) was then combined to determine clus- 
ter membership. That is, both estimates were required 
to be consistent with cluster membership in order for a 
galaxy to be considered a cluster member. 

We used the subset (vast majority) of our spectroscopic 
sample with both optical and near-infrared imaging to 
estimate a) the fraction of spectroscopically-confirmcd 
cluster members excluded by our photometric redshifts 
as a function of morphology and b) the fraction of spec- 
troscopic non-members that are photometric members as 
a function of morphology. Since Pthresh was calibrated 
to include as many cluster members as possible rather 
than to exclude all non-members, the former fraction 
tends to be significantly smaller than the latter. Given 
these fractions, we calculated N m (i), the expected num- 
ber of cluster members of type i that were missed by the 
photometric redshifts; and N c (i), the expected number 
of non-members among galaxies of type i that contami- 
nate the photometric redshift member sample. We then 
used this information to compute the corrected number 



of observed cluster members of type i (iVphotoz (i)), and 
finally, the corrected fraction of cluster members of type 

i (/photoz(«)) : 



AWtoz W = ^obs(i) + N m (i) - N c (i) 



(9) 



-^photoz (toi) = N photoz (E)+N pho t oz (S0)+N photoz (Sp+Irr) 

(10) 

/photoz(i) = Aphotoz(i)/Aphotoz(fo£) (11) 

As in 15.3.11 errors were computed using the iGehrelsl 
(|1986f ) approximation using N p hotoz(«) and N p hotoz(ioi). 
The morphological fractions computed using Equation 
[TT1 are shown as hollow circles in Figure O 

As described in fj2] the area of each cluster with near- 
infrared imaging is somewhat smaller than that imaged 
in the optical. In the classic cosmology, an analysis radius 
of 600 kpc extends beyond the near-infrared imaging for 
cll227-1138. As a result, 19 out of 49 galaxies meeting 
the magnitude and aperture requirements for inclusion in 
the morphological analysis lack near-infrared data. Us- 
ing the WMAP cosmology and an analysis radius of 0.6 x 
i?200> a similar problem occurs for both cll227-1138 (21 
out of 50 galaxies lack near-infrared data) and cll232- 
1250 (14 out of 222 galaxies lack near-infrared data). 
Because the P(z) distributions for photometric redshifts 
computed without near-infrared data are broad, they 
have low P c i U st values, and are somewhat more likely to 
be rejected than galaxies with photometric redshifts com- 
puted using the full filter set. Figure [2] shows that this 
effect does not appear to have systematically skewed the 
morphological fractions for cll227-1138 and cll232-1250 
compared to the other methods employed. 

5.3.3. Statistical Background Subtraction 

In the absence of spectroscopic or photometric red- 
shift information for each galaxy, we can still estimate 
the morphological fractions by statistically correcting the 
number of observed galaxies of a given type to account 
for sources that lie in the field: 



N 



field 



afield 



.4, 



(12) 



(13) 



(14) 



(15) 



iVstat(i) = N ohs (*) - N RcM P(i) 
N sUt (tot) = N obs (tot) - A fic i d 

/stat(i) = Astat(«)/Astat(fo£). 

Here, A is the area of the aperture described in f|57 
Nobs(iot) is the total number of galaxies meeting the 
magnitude criterion ( ^5.1[) within that aperture; N ODS («) 
is the number of these galaxies which have morphology 
i) Nfieid is the number of the observed galaxies that 
are expected to be field members; and N sta t(foi) and 
N s tat(i) are the background-subtracted number of galax- 
ies in the aperture and the number of type i, respec- 
tively. The surface density of field galaxies, Sfi e id, is 
determined by integrati ng the /-band d i fferen tial num- 
ber counts in Table 1 of iPostman et all (|1998f ) down to 
the redshift-dependent /-band apparent magnitude limit 
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adopted for each cluster, as described in §571] We com- 
puted P(i), the fraction of field galaxies of each mor- 
phological type, using data fro m the Medium Deep Sur- 
vey fMDS: lGriffiths et aljl99l . a Hubble Telescope Key 
Project that cataloged the morphologies of intermediate- 
redshift field galaxies down to 7 8 i4 ~ 22 mag. In partic- 
ular, we u sed the classifications of Richard Ellis listed in 
Table 1 of lAbraham et all £1996). The magnitude limits. 
-fauto,iim> of our morphological analysis vary from clus- 
ter to cluster and with cosmology, but generally range 
from 21.3 to 23 mag. The values of P(i) are a function 
of / a uto, linn decreasing to faint magnitudes for early- type 
galaxies and increasing to faint magnitudes for late types. 
For those clusters with J a uto,lim > 22 mag, we adopt the 
values of P(i) down to /gi4 = 22 mag. Based on the 
behavior of P(i) versus Igi4, this procedure likely over- 
estimates P(E), f (SO), and P(Sp) and underestimates 
P(Sp + Irr) and P(Irr). These trends translate into un- 
derestimates of the E, SO, and Sp morphological fractions 
and overestimates of the Sp+Irr and Irr fractions. These 
mis-estimates will be most severe for cll054-1245, cll216- 
1201, and cll354-1230, which have the faintest values of 
^auto,iim- The morphological fractions computed using 
statistical background subtraction and Equation [T5] are 
shown as hollow triangles in Figure [2j Examination of 
this figure does not reveal any obvious biases with re- 
spect to other methods that could be attributable to our 
P(i) estimates. 

As i n §5.3.11 and \ 5.3.2| errors were computed us- 
ing the iGehrelsl (|1986f ) approximation and N sta t(?) and 

N Bt at (tot). 

6. RESULTS 

We have estimated the morphological fractions of 10 
EDisCS clusters at 0.5 < z < 0.8 using four differ- 
ent background-subtraction techniques (absolute limits 
from spectroscopic redshifts, direct estimates from spec- 
troscopic redshifts, photometric redshifts, and statistical 
subtraction), within apertures of different radii (600 kpc 
and 0.6 x P200), an d for two different cosmologies (the 
classic cosmology: Hq = 1, A = 0, Ho = 50 km s _1 
Mpc" 1 ; and the WMAP cosmology: Q = 0.3, A = 0.7, 
H = 70 km s _1 Mpc -1 ). Figure [2] shows how the dif- 
ferent methods compare for an aperture of 0.6 xi? 2 oo 
in the WMAP cosmology. The direct estimate methods 
are consistent with the absolute limits derived from our 
spectroscopy, and with one another. In addition, there is 
no systematic trend for one estimate to produce higher 
or lower fractions than any other. For each cluster we 
therefore adopt a single method: direct estimate from 
spectroscopy, photometric redshifts, or statistical back- 
ground subtraction. The adopted method, which differs 
from cluster to cluster, is chosen in the following way. For 
a given cluster, median E, SO, and Sp+Irr fractions are 
selected from among the estimates of all three methods. 
The method selected most often in this process is the one 
adopted for the cluster. The E, SO, and Sp+Irr fractions 
estimated using the chosen method for each cluster are 
shown in Tables [3] and |4] and are the quantities plotted 
in Figures [3] through [5] 

The morphological fractions in Figure [2] appear to be 
correlated with cluster number, which is determined by 
right ascension. Spearman and Kendall rank correlation 
tests show that marginal correlations are detected at the 



~2<r level for the SO and Sp+Irr fractions. We find no 
correlation between cluster redshift and cluster number 
that could explain this trend. However, we find that the 
cluster velocity dispersion is correlated with the cluster 
number, also at the ~2cr level. Since there is no reason 
for any physical property of the cluster to correlate with 
cluster number, these correlations hint that the morpho- 
logical fractions correlate with cluster velocity dispersion, 
a possibility we revisit in §6.21 

6.1. Evolution of cluster morphological content 

To assess the degree of evolution in the morphologi- 
cal content of rich clusters, we searched the literature 
for similar analyses of clusters spanning a broad redshift 
range. 

The first comparison sample consists of z < 0.5 clusters 
with morphological fractions originally measured by dif- 
ferent groups using a vari ety of methods, but reanalyzed 
in a uniform manner by iFasano et al.l (|2000f ) , hereafter 
F00. The F00 sample includes 55 l ow redshift clusters 
(|Dresslerlll980at iDressler et a!]|1997D: the t en MORPHS 
clusters at 0.37 < z < 0.5 (|Dressler et all 119971 ): three 
clusters at z ~ 0.3 plus A2218 and A1689 at z = 0.18 
(ICouch et al.lll998fk and nine clusters at 0.1 < z < 0.25 
(|Fasano et al.l 120001) . Using this l arge sample o f uni - 
formly analyzed clusters at z < 0.5. lFasano et all (2000) 
confirmed the MORPHs finding that the SO fraction in 
clusters has doubled since z ~ 0.5, while the late-type 
fraction has decreased by a similar factor. 

We also c o mpar e our results with that of 
I Post man et all (|2005l ). who we will refer to as P05. 
They computed the morphological fractions within the 
virial radii of 7 clusters at 0.8 < z < 1.27. This is the 
only large study of clusters at these redshifts in which 
SO galaxies are separately classified. As previously 
discussed, the SO population is an important one to 
track at high redshift. 

In Figure El we show the evolution at z < 0.8 of 
the morphological fractions within rich clusters, using 
EDisCS data in conjunction with the F00 sample de- 
scribed above. All points were computed within 600 kpc, 
using the classic cosmology. We do not include measure- 
ments made on the Postman sample in this plot because 
they were computed within i? 2 oo using the WMAP cos- 
mology. As discussed previously, the F00 sample shows 
a systematic decrease in the SO fraction from 55% at 
z = to 20% at z ~ 0.5. Where the redshift range 
of the EDisCS clusters overlaps with that of F00, the 
morphological fractions derived from the two samples 
are generally consistent. However, the longer redshift 
baseline afforded by the addition of the EDisCS data re- 
veals that the SO fraction is actually flat over the range 
0.4 < z < 0.8. This could mean that z ~ 0.4 is a special 
epoch after which the SO fraction in cluster cores begins 
to grow. The coincidence of this special epoch with the 
redshift where the two samples overlap raises the ques- 
tion of whether the samples were analyzed in different 
ways. However, Section [5] describes our extensive at- 
tempts to control systematics in the EDisCS sample and 
to conform to the analysis presented by F00. Alterna- 
tively, it is possible that we are missing some SO-rich 
clusters at z ~ 0.4 which would reveal that the growth 
of the SO fraction in clusters at z < 0.4 is slower than 
suggested by the present data. Indeed, the scatter in 
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the morphological fractions at z ~ 0.4 is smaller than at 
either higher or lower redshifts. 

The FOO+EDisCS data are consistent with either a 
bend in the morphological fractions at z ~ 0.4 or a 
smooth continuation of the trends observed at z < 0.5, 
but with a flatter overall slope than suggested by the F00 
data alone. More observations at z ~ 0.4 and z > 0.8 
are required to distinguish between these two possibili- 
ties. As discussed above, P05 have analyzed the mor- 
phological fractions, including the SO fraction, in clus- 
ters at z > 0.8. In Figure 0J we show the morphological 
fractions of the EDisCS and P05 clusters as functions of 
redshift. All estimates were made using the WMAP cos- 
mology, but the EDisCS fractions were computed within 
0.6 x i?200 (see N5.2[) . while the P05 fractions were com- 
puted within i?200- P05 find that the E and E+S0 frac- 
tions decrease out to 0.6 x i?200, while the SO fraction 
stays roughly constant and the Sp+Irr fraction increases. 
Between 0.6 x -R200 an d i?200> ah the fractions are flat. 
According to Table 4 of P05, the E, SO, early, and late 
type fractions computed within -R200 are factors of ap- 
proximately 0.79, 0.96, 0.85, and 1.15 times the fractions 
computed within 0.6 x i?200- However, no corrections 
have been made to the points shown in Figure 2J 

Figure [4] shows that the strong evolution in morpho- 
logical fractions seen at z < 0.4 does not continue at 
z > 0.4. Indeed, the E+S0 (Sp+Irr) fraction appears 
to be larger (smaller) in the higher-redshift P05 sample 
compared to the EDisCS clusters. However, a Spearman 
rank correlation analysis shows that there is no statisti- 
cally significant evidence for any evolution over the en- 
tire redshift interval 0.4 < z < 1.25, or in the individual 
EDisCS and P05 samples. Additional clusters in this 
redshift range are necessary to reveal the presence of any 
weak correlation. 

How does the selection requirement for the EDisCS 
clusters to exhibit a red sequence affect the interpreta- 
tion of Figures [3] and Hp If the red sequence takes time 
to build up, it may be expected that this requirement 
selects for clusters that are also dynamically evolved. In 
fact, the EDisCS sample includes clusters that are clearly 
non-spherical, as well as clusters that display significant 
substructure (see H6.2p . The red sequence requirement 
may also be expected to select for clusters with high 
E and/or SO fractions. Figure [3] demonstrates that the 
EDisCS clusters do not contain large E fractions com- 
pared to either the optically-selected F00 sample at lower 
redshifts or the (primarily) X-ray-selected P05 sample at 
higher redshifts. Thus, either such a bias is weak, or it is 
shared by the comparison samples. Given the correlation 
between galaxy color and morphology, any clusters lack- 
ing a red sequence would likely have significantly different 
morphological fractions from the EDisCS, F00, and P05 
samples, leading to increased scatter in Figures [3J and 
01 While both the incidence of such clusters and their 
morphological content are currently impossible to quan- 
tify, in the next subsection we evaluate the dependence 
of morphological fractions on another cluster property: 
velocity dispersion. 

6.2. Correlation between Morphological Fractions and 
Cluster Velocity Dispersion 

In the previous subsection, we argued that there is 
no systematic evolution in the morphological fractions 



within rich clusters at 0.4 < z < 1.25. Figures [3] and [4] 
indicate that the scatter in these fractions is large. How 
much of this scatter is due to a correlation between mor- 
phological fractions and cluster mass? Such a correlation 
may be expected in either nature or nurture scenarios. 
In the former it is due to the fact that the most massive 
clusters collapsed at earlier times. In the latter it could 
be due to a higher efficiency of morphological transfor- 
mations in more massive clusters. 

P05 found that the E, SO, and E+S0 fractions within 
-R200 of seven z ~ 1 clusters increase with increasing 
bolometric cluster X-ray luminosity, although the corre- 
lations are significant only at the <3cr level. However, 
they find no correlation between the E+S0 fraction and 
the X-ray temperature or the cluster velocity dispersion, 
perhaps because of small number statistics. In Figure [5] 
we plot the morphological fractions of the EDisCS and 
P05 samples against cluster velocity dispersion. From 
this figure it is apparent that clusters with larger velocity 
dispersions harbor a higher fraction of early type galax- 
ies and fewer late type galaxies. Spearman and Kendall 
rank correlation tests show that these visual impressions 
are statistically significant at the <3cr level. 

While it is tempting to interpret this correlation as one 
between morphological content and cluster mass, the ve- 
locity dispersion of a cluster is directly related to its mass 
only if the cluster is virialized. The velocity dispersion 
may overpredict mass if the cluster is experiencing signif- 
icant merging events. The degree of substructure in five 
of the clusters in this work (cll040-1155, cll054-1146, 
cll054-1245, cll21 6-1201, and cll232-1 250) has been 
studied in detail by Hallidav et alj (12004D usi ng Dressler- 
Shectman tests ([Dressier fc Shectmanl fl988l) . They de- 
tect substructure in cll232-1250 and cll216-1201 with 
more than 95% confidence. For cll040-1155 and cll054- 
1245, not enough spectra were available to provide firm 
evidence for substructure. No evidence for substructure 
was found for cll054-1146. They note that the two clus- 
ters showing clear evidence of substructure also have the 
largest velocity dispersions in the EDisCS sample, and 
caution against using the velocity dispersions for these 
two clusters as a proxy for mass. Further analysis on the 
remaining EDisCS clusters is ongoing (Milvang-Jensen 
et al., in preparation). 

What is the situation for the z < 0.5 sample, 
where evolution is observed? Velocity dispersions 
were available in th e literature for 14 of the clusters 
in the F00 sample (ICouch fc Sharpies! Il987t iGudehusI 
19891 iGirardi fc Mezzettl FoOl: De Propris et alj 120021 : 
Bettoni et al.l 120061 ). Figure [S] shows the morphological 
fractions for these 14 F00 clusters as a function of clus- 
ter velocity dispersion. Without a larger number of data 
points it is difficult to make a conclusive statement re- 
garding the existence of a correlation. Spearman and 
Kendall rank correlation tests indicate that none of the 
morphological fractions are significantly correlated with 
cluster velocity dispersion. However, removal of two clus- 
ters (A389 and A3330) with unusually high SO fractions 
for their velocity dispersions results in a 2a detection 
of a correlation between SO fraction and cluster velocity 
dispersion. Additional velocity dispersions for the F00 
sample would greatly aid an assessment of any correla- 
tion, which is necessary for understanding the observed 
trends between morphology and redshift. 
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Comparing the x-axes of Figures [5] and [6l we see that 
the F00 sample at z < 0.5 includes some clusters with 
very high velocity dispersions (a > 1200 km s^ 1 ), while 
the z > 0.5 EDisCS+P05 sample does not. Figure [7] 
shows the cluster velocity dispersions of the EDisCS, 
F00, and P05 samples as a function of lookback time. 
The clusters with the highest velocity dispersions in these 
samples lie at low redshift. Unfortunately, velocity dis- 
persions are unavailable for a significant fraction of the 
low redshift sample. These are urgently needed to under- 
stand the extent to which sample selection is responsible 
for the apparent evolution in the morphological content 
of clusters. 

7. DISCUSSION AND CONCLUSIONS 

Using high resolution imaging afforded by the ACS in- 
strument on board the Hubble Space Telescope, we have 
shown that the morphological content of rich clusters 
with velocity dispersions in the range a = 200 — 1200 
km s^ 1 shows no evidence for evolution over the redshift 
range 0.4 < z < 1.25 (see Figures [3] and 0} . Although 
the scatter is significant, typical morphological fractions 
for clusters in this redshift range are 0.3, 0.15, and 0.55 
for E, SO, and Sp+Irr galaxies, respectively. In contrast, 
studies of clusters at lower redshifts have shown that the 
elliptical fraction remains constant between z = and 
z ~ 0.5, while the SO fraction decreases by a factor of ^2 
and the Sp+Irr fraction increases by a similar amount 
over the same redshift range (e.g. F00). Our data show 
that the observed evolution in cluster SO and Sp+Irr pop- 
ulations does not continue beyond z ~ 0.4, at least for 
the velocity dispersions probed in this study. 

How do our results concerning the global evolution of 
galaxy morphologies within rich clusters translate into 
a statement regarding the evolution of the morphology- 
density relation? Although the centers of cluster cores 
are regions of high local galaxy density, our morphologi- 
cal fractions were computed within sizable apertures (see 
H5.2[) . so the average environment we are probing is of 
moderate density. For example, the average galaxy sur- 
face density within a radius of 0.6 x i?2oo (^o = 0.3, 
A = 0.7, H = 70 km s" 1 Mpc" 1 ) for the EDisCS clus- 
ters used in this analysis ranges from ~40-175 Mpc -2 . 
The morphological fractions of the F00 and P05 clusters 
were likely computed using galaxies inhabiting environ- 
ments of similar average density. These results may point 
to a lack of evolution in the morphology density relation 
between z = 0.8 and z = 0.4, with a subsequent in- 
crease in the SO population and a decrease in the Sp+Irr 
population. Larger samples in this redshift range are 
needed to understand the scatter in the morphological 
fractions, and to rule out weak evolution. However, our 
current data are consistent with rec ent studies of the 
morphology-density relation at z ~ 1 (jSmith et al.ll2005t 
iPostman et~alll2005| ). 

Several studies of nearby galaxies indicate that the ob- 
served relation between morphology and environment is 
a reflection of a more primary relation between star for- 
mation rate (SFR) and environm ent (Kauff mann et all 
[200l iChristlein fc Zabludoffl [20051 : [Blanton et al.ll2005tT 
The star formation properties of cluster, group, and field 
galaxies in the EDisCS spectroscopic sample have been 
measured and compared to a lo w redshift sample from 
the Sloan Digital Sky Survey by iPoggianti et all (|2006D 



(P06). Based on these data, P06 put forth a model in 
which passive galaxies devoid of ongoing star formation 
are made up of two separate populations: "primordial" 
galaxies whose stars all formed at z > 2.5 and galax- 
ies whose star formation lasted until later times but was 
ultimately "quenched" due to entering a cluster-like en- 
vironment (> 1O 14 M0). In this model, the primordial 
galaxies may be identified with the elliptical galaxies that 
make up ^30% on average of the galaxy populations in 
clusters at z < 1, and perhaps some of the SO galaxies 
(see below). This identification is consistent with both 
the lack of evolution in the elliptical fraction at z < 1 
and with the old ages inferred for elliptical stellar popu- 
lations. It is tempting to identify the quenched galaxies 
with SOs. However, the fraction of quenched galaxies in 
the model increases at z < 0.8, while the SO fraction 
only increases at z < 0.4. These timescales are consis- 
tent if it takes roughly a billion years for a galaxy to 
resemble an S O after the cessation of star formation in a 
spiral galaxy (iDressler et aT. 1999; P oggianti et al.lfl999l : 
iTran et al.ll2003f h However, the formation of the SO pop- 
ulation observed in clusters at 0.4 < z < 1 cannot be 
explained in this way. Perhaps these are "primordial" 
SOs that assumed an SO morphology upon or soon after 
formation at z > 2.5. If so, the SO galaxies in clusters at 
z > 0.4 should be, on average, older and more massive 
than cluster SOs at z < 0.4. Another possibility is that at 
least some late types may transform into SOs in systems 
less massive than 10 14 M Q , perhaps by an entirely dif- 
ferent mechanism than operates in > 10 14 M(T) systems . 
The existence of SOs in groups (e.g. iHickson et al.|[l989l ) 
suggests that this may be the case. 

If a significant fraction of SO galaxies in nearby clus- 
ters were quenched, what did they look like beforehand? 
The parallel decline in the Sp+Irr fraction as the SO frac- 
tion increases between z = 0.4 and z = suggests that 
some subset of late type galaxies transformed into SOs 
subsequent to quenching. For the EDisCS clusters, we 
computed the fractions of Sp and Irr galaxies separately 
(see Tables [3] through 0]) . As with the other morpholog- 
ical fractions at z > 0.4, the Sp and Irr fractions do not 
vary systematically with redshift. Moreover, the Irr frac- 
tions in the EDisCS clusters are very small. We conclude 
that Sp galaxies and not Irr galaxies are the precursors 
of the quenched SOs in local clusters. 

In the P06 model, both nature and nurture (within 
clusters) play a role in the environmental dependence 
of star formation and therefore morphology. The nur- 
ture component within clusters is consistent with the ob- 
servation that the relation between star formation and 
environment in nearby clusters is not solely a reflec- 
tion of the relation between stellar age and environ- 
ment, as would be expected in a pure nature scenario 
(|Christlein fc Zabludoffl [2005h . The P06 model is also 
compatible with observations which show that the SFR- 
density relation out to z ~ 1 extends to very low local 
densities, comparable to those found at the virial radius 
of clusters and even outside clusters dLewis et al.l 120021 : 
iGomez et all2003l: ICooper et al.ll2006D . The low density 
relations could be set up by the dependence of galaxy 
properties on initial conditions, or perhaps a different 
mechanism creates SO galaxies in low-density environ- 
ments, or a combination. 

Although the model forwarded by P06 is an attrac- 
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tive framework within which to view our observational 
results on the morphological content of rich clusters, it 
is not necessarily the only model which could explain 
all of the data accumulating on the evolution of galaxy 
properties with environment. Furthermore, it does not 
identify a mechanism responsible for the "quenching" , 
although it does specify a timescale (3 Gyr) and a mass 
scale (10 14 Mq) of a workable mechanism. 

In addition to analyzing the evolution of the morpho- 
logical content in galaxy clusters out to z = 1.25, we have 
also studied the dependence of the morphological con- 
tent on the cluster velocity dispersion (see Figures [5] and 
[6]). We find that the early and late- type fractions in the 
EDisCS and P05 cluster samples correlate with cluster 
velocity dispersion at a statistically significant level. This 
correlation highlights the importance of understanding 
global cluster properties in samples where evolution is 
observed, such as the F00 sample. Unfortunately, only 
limited velocity dispersion information is available for the 
F00 clusters. The existing information for 14 clusters in- 
dicates that they have a higher average velocity disper- 
sion than either the EDisCS or P05 clusters. In addition, 
the morphological fractions of these 14 clusters do not 
correlate strongly with cluster velocity dispersion. It is 
possible that additional velocity dispersions for the F00 
sample will reveal an underlying trend, that additional 
studies at high redshift will prove the trend observed 
in the higher redshift EDisCS+P05 sample to be spuri- 
ous, or that the relation between morphological fraction 



Abadi, M. G., Moore, B., & Bower, R. G. 1999, MNRAS, 308, 947 
Abraham, R. G., Tanvir, N. R., Santiago, B. X., Ellis, R. S., 

Glazebrook, K., & van den Bergh, S. 1996, MNRAS, 279, L47 
Bckki, K., Couch, W. J., & Shioya, Y. 2002, ApJ, 577, 651 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Bettoni, D., Moles, M., Kjasrgaard, P., Fasano, G., & Varela, J. 

2006, A&A, 452, 811 
Blanton, M. R. et al. 2001, AJ, 121, 2358 

Blanton, M. R., Eisenstein, D., Hogg, D. W., Schlegel, D. J., & 

Brinkmann, J. 2005, ApJ, 629, 143 
Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 476 
Byrd, G., & Valtonen, M. 1990, ApJ, 350, 89 
Christlein, D., & Zabludoff, A. I. 2005, ApJ, 621, 201 
Clowe, D. et al. 2006, A&A, 451, 395 
Cooper, M. C. et al. 2006, MNRAS, 370, 198 

Couch, W. J., Barger, A. J., Smail, I., Ellis, R. S., & Sharpies, 

R. M. 1998, ApJ, 497, 188 
Couch, W. J., & Sharpies, R. M. 1987, MNRAS, 229, 423 
Dalcanton, J. J. 1996, ApJ, 466, 92 
De Propris, R. et al. 2002, MNRAS, 329, 87 
Dressier, A. 1980a, ApJs, 42, 565 
— . 1980b, ApJ, 236, 351 
Dressier, A. et al. 1997, ApJ, 490, 577 
Dressier, A., & Shectman, S. A. 1988, AJ, 95, 985 
Dressier, A., Smail, I., Poggianti, B. M., Butcher, H., Couch, W. J., 

Ellis, R. S., & Oemler, A. J. 1999, ApJS, 122, 51 
Fasano, G., Poggianti, B. M., Couch, W. J., Bettoni, D., 

Kjasrgaard, P., & Moles, M. 2000, ApJ, 542, 673 
Finn, R. A. et al. 2005, ApJ, 630, 206 
Gomez, P. L. et al. 2003, ApJ, 584, 210 
Gehrels, N. 1986, ApJ, 303, 336 
Girardi, M., & Mezzetti, M. 2001, ApJ, 548, 79 
Gonzalez, A. H., Zaritsky, D., Dalcanton, J. J., & Nelson, A. 2001, 

ApJS, 137, 117 
Goto, T. et al. 2002, PASJ, 54, 515 

Goto, T., Yagi, M., Tanaka, M., & Okamura, S. 2004, MNRAS, 
348, 515 

Griffiths, R. E. et al. 1994, ApJ, 435, L19 

Gudehus, D. H. 1989, ApJ, 340, 661 

Gunn, J. E., & Gott, J. R. I. 1972, ApJ, 176, 1 



and cluster velocity dispersion is itself a function of red- 
shift. Larger cluster samples with robust velocity disper- 
sions and morphologies are needed to determine which 
of these options is most likely. If the morphology-density 
relation is driven by environmental processes in clusters, 
such measurements are essential for determining how the 
efficiency of these transformations depends on cluster ve- 
locity dispersion, which, as discussed in ^6.21 is related to 
both the cluster mass and its dynamical state. 



We greatly appreciate the timely and cheerful assis- 
tance of Galina Soutchkova in planning the HST obser- 
vations. We also thank the referee, Manolis Plionis, for 
feedback which improved the paper. VD acknowledges 
funding from the Graduate Student Researchers Pro- 
gram. This work was supported by NASA grant HST- 
GO-09476.01. JJD was partially supported by the Alfred 
P. Sloan Foundation. The Medium Deep Survey catalog 
is based on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. The Medium-Deep Survey anal- 
ysis was funded by the HST WFPC2 Team and STScI 
grants G02684, G06951, G07536, and G08384 to Prof. 
Richard Griffiths and Dr Kavan Ratnatunga at Carnegie 
Mellon University. 



Halliday, C. et al. 2004, A&A, 427, 397 

Hickson, P., Kindl, E., & Auman, J. R. 1989, ApJS, 70, 687 

Icke, V. 1985, A&A, 144, 115 

Kauffmann, G., White, S. D. M., Heckman, T. M., Menard, B., 
Brinchmann, J., Chariot, S., Tremonti, C, & Brinkmann, J. 

2004, MNRAS, 353, 713 

Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 
692 

Lavery, R. J., & Henry, J. P. 1988, ApJ, 330, 596 

Lewis, I. et al. 2002, MNRAS, 334, 673 

Melnick, J., & Sargent, W. L. W. 1977, ApJ, 215, 401 

Mihos, J. C. 2004, in Clusters of Galaxies: Probes of Cosmological 

Structure and Galaxy Evolution, 277 — h 
Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 
Poggianti, B. M. 1997, A&AS, 122, 399 

Poggianti, B. M., Smail, I., Dressier, A., Couch, W. J., Barger, 
A. J., Butcher, H., Ellis, R. S., & Oemler, A. J. 1999, ApJ, 518, 
576 

Poggianti, B. M. et al. 2006, ApJ, 642, 188 
Postman, M. et al. 2005, ApJ, 623, 721 

Postman, M., Lauer, T. R., Szapudi, I., & Oegerle, W. 1998, ApJ, 
506, 33 

Quilis, V., Moore, B., & Bower, R. 2000, Science, 288, 1617 
Richstone, D. O. 1976, ApJ, 204, 642 
Rudnick, G. et al. 2001, AJ, 122, 2205 
— . 2003, ApJ, 599, 847 

Smail, I., Dressier, A., Couch, W. J., Ellis, R. S., Oemler, A. J., 

Butcher, H., & Sharpies, R. M. 1997, ApJS, 110, 213 
Smith, G. P., Treu, T., Ellis, R. S., Moran, S. M., & Dressier, A. 

2005, ApJ, 620, 78 

Thomas, T., & Katgert, P. 2006, A&A, 446, 31 
Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 
Tran, K. H., Simard, L., Illingworth, G., & Franx, M. 2003, ApJ, 
590, 238 

Treu, T., Ellis, R. S., Kneib, J., Dressier, A., Smail, I., Czoske, O., 

Oemler, A., & Natarajan, P. 2003, ApJ, 591, 53 
White, S. D. M. et al. 2005, A&A, 444, 365 

Zaritsky, D., Nelson, A. E., Dalcanton, J. J., & Gonzalez, A. H. 

1997, ApJ, 480, L91 
Zucca, E. et al. 2006, A&A, 455, 879 



EDisCS Morphological Content 



11 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 

R R 






0.4 r 




0.3 ( 


>< 






0.2 \ 








0.1 ■ 










0.0 0.1 0.2 0.3 0.4 0.5 

t 



0.0 0.1 0.2 0.3 0.4 0.5 

X 2 



Fig. 1. — Comparison of parameters derived from ID profile fits to visually-classified ellipticals and SOs down to Tauto = 23 mag in 
cll216-120f. From top to bottom, this figure shows the fractional differential (left) and cumulative (right) distribution functions of the 
width of the ellipticity distribution within a given galaxy (R), the bulge-to-total fraction (B/T), and the \ 2 goodncss-of-fit parameter 
for a pure de Vaucouleurs component (see [2}. The solid lines indicate the distributions for ellipticals, and the dotted lines indicate the 
distributions for SO galaxies. The significance levels of the two-sided Kolmogorov-Smirnov tests comparing the E and SO distributions are 
shown in the bottom right of the cumulative distribution plots. The R and B/T distributions for Es and SOs are different at very high 
significance levels, while the \ 2 goodness-of-fit parameters are different at a low significance level. These results show that galaxies that 
are classified as Es are structurally distinct from those classified as SOs. 
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Fig. 2. — A comparison of four methods for estimating the E, SO, Sp, and Sp^ 
shaded grey regions span the full range allowed by our spectr oscopic data, while the squares show the best estimates of the morphological 
fractions determined using spectroscopic redshifts (see £|5. 3.11 1. Circles and triangles represent the morphological fractions determined using 
photometric redshifts (see £|5. 3.21 1 and statistical background subtraction (see i|5. 3.31 1. respectively. All calculations were performed using 
the WMAP cosmology and an aperture of 0.6 * J?200- 
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Fig. 4. — The evolution of the E, SO, E+S0, and Sp+Irr fractions as traced by EDisCS clusters (filled circles) and P05 clusters (hollow 
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Fig. 5. — Morphological fractions as a function of cluster velocity dispersion for the EDisCS clusters (filled circles) and the subset of the 
P05 cluster sample for which velocity dispersions are available (hollow circles). The EDisCS fractions were computed within a radius of 
0.6 X -R200 and the P05 fractions were computed within a radius of i?200- All fractions were computed using the WMAP cosmology. 
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measurements for ten of these clusters come from Girardi & Mczzctti (2001). The remainder come from [Bettoni et' al. (2006) (A3330 
with a = 695 km s" 1 ), IDe Propris et al.l (I2002T) (A389 with a = 667 km s" 1 ), IGudehusi (ITMSI) (A1689 with a = 1800 km s _1 ), and 
ICouch fc Sharpies (1987) (AC118 with a = 1950 km s — 1 ). Both A3330 and A389 have large SO fractions compared to other clusters with 
similar velocity dispersions. In addition, A1689 and AC118 are the F00 clusters with the largest velocity dispersions, and largely drive 
the correlation between morphological fractions and velocity dispersion. This plot illustrates the need for homogeneous studies of velocity 
dispersions for clusters where evolution in morphological fractions is observed. 
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TABLE 1 
The EDisCS HST Sample 



Number 


Name 


RA (J2000) 


Dec (J2000) 
(»:':") 


Z <T -^200,1 -f?200,2 






(h:m:s) 


(km s" 1 ) (Mpc) (Mpc) 



1 


cll037-1243 


10 


:37:51.4 


-12; 


:43: 


:26. 


.6 


0. 


58 


319±g 





.56 





.57 


2 


C11040-1155 


10 


:40:40.3 


-11: 


:56 


:04. 


.2 


0. 


70 







.65 





.70 


3 


cll054-1146 


10 


: 54:24.4 


-11; 


:46: 


:19. 


.4 


0. 


70 


589+to 





.92 





.99 


1 


C11054-1245 


10 


: 54:43. 5 


-12: 


:45: 


:51. 


.9 


0. 


75 


504t^ 3 





.75 





.82 


5 


clll03-1245b 


11 


:03:36.5 


-12; 


:41: 


:22 


.3 


0. 


70 


252ti 





.39 





.42 


6 


C11138-1133 


11 


:38:10.2 


-11: 


:33: 


:37. 


.9 


0. 


.48 


732^ 


1 


.41 


1 


.40 


7 


cll216-1201 


12 


:16:45.3 


-12: 


:01: 


:17. 


.6 


0. 


79 


1018+ 77 


1 


.47 


1 


.61 


8 


C11227-1138 


12 


:27:58.9 


-11: 


:35: 


:13. 


.5 


0. 


64 







.95 


1 


.00 


9 


cll232-1250 


12 


:32:30.3 


-12: 


:50: 


36. 


.4 


0. 


54 


losol^ 9 


1 


.95 


1 


.99 


10 


cl!354-1230 


13 


: 54:09. 8 


-12: 


:31: 


:01. 


.5 


0. 


76 







.96 


1 


.05 



Note. — The RA, Dec coordinates are the positions of the brightest clust er galaxy. The 
redshi fts (2) and line-of-sight cluster velocity dispersions (cr) were taken from Halliday ct al. 
(2004) and Milvang- Jensen et al. (in preparation). Th e values of the virial ra dius, -R20O1 listed 
in columns 6 and 7 were computed using Equation 1 in Poggianti ct al. (2006) using the classic 
cosmology and the WMAP cosmology, respectively, as described in [[T] 



TABLE 2 

Notes on Parameters in Morphological Catalogs. 



Heading 



Description 



ID 

RA 
DEC 

^auto 

Type 



SO disk flag 

bar flag 
edge-on flag 
small flag 
LSB flag 
defect flag 

dust flag 
disturbance flag 
comments 



EDisCS ID 

Right Ascension in decimal degrees (J2000) 
Declination in decimal degrees (J2000) 

Total / magnitude; SExtractor MAG-AUTO parameter measured from ground-based /-band 
images 

star=-7, non-stellar but too compact to see structure=-6, E=-5, S0=-2, Sa=l, Sb=3, 
Sc=5, Sd=7, Sm=9, Irr=ll, no HST data corresponding to ground-based object=lll, 
unclassifiable=66 

If any of the classifiers noted that the B/D ratio is small but the disk is featureless, this flag 
is 1. 

If any of the classifiers noted the presence of a bar, this flag is 1. 

If any of the classifiers noted that this galaxy is edge-on, this flag is 1. 

If any of the classifiers noted that this galaxy is small, this flag is 1. 

If any of the classifiers noted that this galaxy is low-surface-brightness, this flag is 1. 

If any of the classifiers noted that the image of this galaxy was defective, due for example to 

cosmic rays or incomplete coverage, this flag is 1. 

If any of the classifiers noted the presence of dust in the galaxy, this flag is 1. 
If any of the classifiers noted that this galaxy is disturbed, this flag is 1. 
Additional comments by any of the classifiers 



TABLE 3 

Morphological Fractions within EDisCS Clusters: R = 600 kpc, fio = 1, A = 0, Ho = 50 km 

s _1 Mpc -1 



Cluster Name 


fa 


/so 


/e+so 


/Sp+lrr 


fSp 


fin 



cll037-1243 


0. 


.2811° 


124 
.140 





OOOlo 


109 
.000 


0. 


.281 


+0 
-0 


124 
146 


0. 


.719 


+0. 
-0 


146 
124 





.625+° 


.138 
.156 


0. 


.095+° 


174 

.050 


cll040-1155 


0. 


.377 +0 


130 
.116 





066lg 


003 
.058 


0. 


.444 


+0. 
-0 


141 
123 


0. 


.556 


+0. 
-0 


128 
141 


0. 


.419+° 


141 
.115 


0. 


1371° 


127 
.071 


cll054-1146 


0. 


.2451° 


071 
.060 





ooolg 


031) 



0. 


.245 


+0 
-0 


071 
060 


0. 


.755 


+0 
-0 


1160 
071 





.7551° 


060 
.071 





OOOIq 


036 
.000 


cll054-1245 





.300 + o 


107 

.090 





267lg 


104 

.087 


0. 


.567 


+0. 
-0 


102 
108 


0. 


433 


+0. 
-0 


108 
102 





.433+° 


108 
.102 





OOOIq 


1)01) 
.000 


ell 103- 1245b 


0. 


.250+° 


120 

-OKI) 





OOO+o 


070 
.000 


0. 


250 


+0 
-0 


120 
080 


0. 


.750 


+0 
-0 


oso 
120 





.7501° 


DM) 
.120 





OOOIq 


070 
.000 


C11138-1133 





■ 3051° 


164 

.120 





■0951° 


113 

0S4 


0. 


.400 


+0. 
-0 


162 
145 


0. 


.600 


+0. 
-0 


145 
162 





600l° 


145 

.162 





OOOIq 


115 

.000 


C11216-1201 


0. 


490lo 


030 
.020 





.220+0 


1121) 
.020 


0. 


.710 


+0 
-0 


020 
020 





.290 


+0 
-0 


020 
020 





.270+°, 


1)21) 
.020 





.020+0 


010 
.010 


C11227-1138 





.290+0 


105 
.136 





.1461° 


.157 
.005 


0. 


.436 


+0. 
-0 


160 
162 


0. 


.564 


+0. 
-0 


162 
160 





3941° 


167 
.164 





■ 1701° 


132 
.120 


C11232-1250 


0. 


.350+° 


040 
.040 





.1701° 


030 
.0.80 





.530 


+0 
-0 


10 
040 


0. 


.470 


+0 
-0 


IO 
040 





.470+0 


oil) 

.040 





OOOIq 


010 



cll354-1230 





■ 1701° 


070 

.050 





290lg 


070 

.060 


0. 


.450 
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Note. — The morphological fractions for a given cluster were computed using the background 
subtraction method adopted for that cluster (see J6}. 
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TABLE 4 

Morphological Fractions within EDisCS Clusters: R = 0.6 x -R200, Oo = 0.3, A = 0.7, 

H = 70 km s _1 Mpc -1 
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Note. — The morphological fractions for a given cluster were computed using the background 
subtraction method adopted for that cluster (see 



